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Abstract 

Using a simulation technique introduced recently, we study winding 
clusters in percolation on the torus and the Mobius strip for different 
aspect ratios. The asynchronous parallelization of the simulation makes 
very large system and sample sizes possible. Our high accuracy results are 
fully consistent with predictions from conformal field theory. The numer- 
ical results for the Mobius strip and the number distribution of winding 
clusters on the torus await theoretical explanation. To our knowledge, 
this study is the first of its kind. 

1 Introduction 

Since the seminal article of Langlands et al. [Q, percolation has enjoyed a 
renaissance in the last decade. Cardy's result |2] for the probability of a crossing 
cluster on a rectangle was the first major breakthrough of conformal field theory 
in critical percolation. Related to this work, Pinson calculated the probability 
of clusters of particular winding numbers in [3] , using the partition sum for the 
torus formulated earlier by di Francesco et al. 4 . 

In 1996 Hu and Lin showed numerically [5] that the probability of more 
than one percolating cluster is non-vanishing, which was proven by Aizenman 
|H] shortly afterwards. In fact, Cardy was able to calculate the asymptotic prob- 
ability of n distinct, simultaneously crossing clusters [x] exactly. One expects 
similar behavior for the number of distinct, simultaneously winding clusters on 
the torus. 

In this article we present numerical results for winding clusters on the torus 
and on the Mobius strip with different aspect ratios, providing strong numerical 
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support for the claim of conformal invariance in two dimensional percolation. 
Apart from the numerical results for r = 1 by Langlands et al. presented in [3], 
the only other numerical studies of this kind, to our knowledge, is one by Ziff 
et al. |5] , who measure the probability of a cross topology for very small system 
sizes (16 x 16) and different twists on a torus with aspect ratio 1, and those 
by Newman and Ziff |SJ ^| , who measure the probability of winding clusters 
on a torus with aspect ratio 1 and system sizes up to 256 x 256 and varying 
occupation probability. Also, for the Mobius strip we are only aware of studies 
on the Ising model [111 I12| . 

1.1 Critical percolation on the torus 

We treat site and bond percolation as the two limiting cases of the more general 
site-bond percolation: Sites are occupied with probability p <s) and bonds are 
activated with probability p (b) . In site percolation all bonds are active, while 
in bond percolation all sites are occupied. Two sites are connected if the bond 
between them is active and both sites are occupied. Two sites belong to the same 
cluster if they are connected by a path along occupied sites and active bonds. 
On a torus and the Mobius strip, these paths may wind around the lattice. 
By our convention, (a, b) counts the number of windings in the horizontal and 
vertical directions, respectively. Vertical and horizontal directions are fixed 
by the definition of aspect ratio as the vertical circumference (waist) over the 
horizontal circumference (ignoring the distortion), see Fig. ^ 

On the torus, it is a topological fact that if a = (b = 0) then the only 
path which is not homotopic to a point and does not intersect itself has b = 1 
(a = 1). There is no topological difference between (a, 6) and (—a, —b) — the 
direction in which these paths are taken is simply inverted. If a ^ and 6^0 
then a and 6 must be relative prime if the path does not intersect itself. Fig. [21 
shows an example for (2, 3), other examples can be found in 0j. 

In the following, winding numbers (a, b) will be called "normalized" , if a = 1 
for b = 0, b = 1 for a = 0, a and b are relative prime if a =/= and b =/= 0, and 
a > 0. 

In addition to the above, we also mention the path with the so-called "cross 
topology" @|. Such a path is produced by intersecting a (1,0) and a (0, 1) path. 
The seam shown in Fig. [I]as a dashed line forms a path with a cross topology. 

To transfer the notion of winding from paths to clusters, one considers all 
(topologically different) paths within a cluster. If there is a path with non-zero 
winding numbers that crosses only paths which are homotopic to a point or 
have the same winding numbers, the cluster is assigned the winding numbers of 
the path. If there is no winding path at all, then this cluster itself is said to be 
homotopic to a point. It transpires that in any other case the cluster contains 
a cross-topological path and the cluster is then said to have a cross topology 
itself. 

The universal probability to obtain n clusters with winding numbers (a, b) 
at aspect ratio r is denoted in the following by V((a. b), n, r). Pinson has derived 
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Figure 1: A torus with r < 1. The winding directions of (a, b) are marked as 
such and the seam is indicated as dashed line. Below, the "unwrapped" lattice 
shows the definition of the aspect ratio, r = l v /lh- 



the formula [3] 



where 



and 



V((a,b),>l,r)=Y / Za3i,b3i (\;r) - (1) 



iez ^ ' zez '* 

iez ^ ' ;ez v '* 



Z m , n (g;r) = ^ 2(e _ 2 , r) exp ^-7rg(— + n 2 r) j , (2) 

OO 

P((o,6),>n,r) = X)P(M),i,r) . (3) 
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Figure 2: A (2,3) winding on a torus with aspect ratio r < 1. The seam is 
shown as a dashed line. 



Here 77(g) is the Dedekind eta function 

oo 

v{q)=q 1,2A \{{l-q k ) ■ (4) 

fe=l 

Correspondingly, the probability of a cross topology is denoted by P(X, r), 
and the probability that all clusters in a given configuration are homotopic to 
a point is denoted by V(0, r). The exact expressions for these probabilities are 

V(X, r) = V(0, r) = i(zc(|,l;r)-Z c (|,i;r)) (5) 

with 

z c {g,f;r) = f ^2 z fm,fn(g;r) . (6) 

m,n£Z 

In the following, we will distinguish exact results from numerical results by 
putting a hat on all exact quantities. Where necessary, the numerical results 
will also have an index indicating the system size and a superscript (a) for site 
percolation and <b> for bond percolation. For example, 'P^^ =3000 2 ((1, 2), 1, 9) 
is the fraction of cases in which we observed a single (1,2) cluster in bond 
percolation with aspect ratio 9 and 3000 2 sites. 

The topological considerations in this paper are mainly technically motivated 
and rather heuristic. For rigorous proofs, we refer to the standard literature 

Multiple, distinct clusters with the same winding number can coexist with- 
out intersecting [I]. This, however, does not apply to a cluster with a cross 
topology: there can only be one such cluster on a torus. As explained below in 
detail, winding clusters with incommensurable winding numbers cannot coexist. 
Thus an entire configuration of the lattice on the torus is characterized by the 
winding numbers of the winding clusters, if any, and their total number. 
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The Mobius strip is a little more complicated in this respect. First of all, 
since the Mobius strip is a non-orientable surface, a reasonable definition of a 
spanning cluster connecting one side to the other is not possible, unlike, for 
example, clusters on a cylinder connecting top and bottom ^^[7J- Winding 
clusters behave rather surprisingly. A single winding cluster, winding around 
only once, is possible. Meanwhile, if two winding clusters coexist, then at least 
one of them must have winding number 2. In general, n winding clusters require 
at least n — 1 clusters with winding number 2. In fact, all winding clusters are 
fully defined by the total number a of windings alone: If a is even, there are 
a/2 winding clusters with winding number 2. If a is odd, there are (a — l)/2 
winding clusters with winding number 2, and 1 cluster with winding number 1. 
Fig. |21 shows a Mobius strip with 2 winding paths adding up to a total winding- 
number 3. The aspect ratio of the Mobius strip is defined as the width of the 
band over the path length along the band, as shown in Fig. □ 



I 




w 



I 

Figure 3: A Mobius strip with two non-intersecting winding clusters, one with 
winding number 2 (full line) and one with winding number 1 (dashed line). 
Below, the "unwrapped" lattice shows the definition of the aspect ratio, r = w/l. 



One might be tempted to extend the study of winding numbers to the Klein 
bottle. In this case, however, the wrapping around the waist is ill-defined, 
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just like for the spanning cluster on the Mobius strip discussed above. The 
twist in the boundary conditions breaks translational invariance; for example, 
a winding path, for which the algorithm (see below) would indicate winding 
numbers (1,0), can be smoothly transformed into a path with winding numbers 
(1, 1). This also makes the Klein bottle fundamentally different from the torus 
with a twist, which has a fixed "offset" for its winding numbers. Therefore, we 
did not study the Klein bottle. 

2 Method 

The simulation is based on a method described in detail in JS] , which in princi- 
ple relaxes all the standard constraints in numerical simulations of percolation, 
such as CPU power, memory, and network capacity. The method is especially 
suited for calculating cluster size distributions and crossing probabilities. Here 
we have used the method to calculate probabilities of winding clusters on the 
torus and the Mobius strip. 

The algorithm runs asynchronously across a cluster of computers (slaves), 
which supply a central node (master) with "patches". These patches are pro- 
duced by simulating percolation on squares and encoding in the boundary in- 
formation about all clusters touching the boundary. The data representation 
is due to Hoshen and Kopelman [FS\ (HK), and the reduction of the data to 
the boundary is essentially a form of "Nakanishi label recycling" (171 ITS] . The 
patches are "glued together" by the master node to form larger lattices with 
various boundary conditions, such as periodic boundaries used here. One of the 
key features of the algorithm is that the same patches can be used for different 
topologies (such as the torus and the Mobius strip) and aspect ratios. To re- 
duce correlations between different results, patches can be randomly permuted, 
mirrored and rotated. We have applied this technique for each aspect ratio. 

Technical details concerning the method can be found in [15| and, especially 
regarding the detection of wrapping clusters, in 14 j. Here we assume that 
a master node has received a number of square patches of size L 2 , and has 
(randomly) glued them together to form a lattice with open boundaries of size N 
and aspect ratio r. At this stage, sites are still encoded in HK style: Every site on 
the border has a pointer attached to it, which points either to its representative 
or to itself. If it points to itself, the site is called a root site. A representative 
may or may not point to another representative, so that the pointers form a tree- 
like structure. The algorithm ensures that a tree has exactly one root and that 
all sites within the same cluster belong to the same tree. Thus, all information 
about a cluster can be stored at the root site and is accessible from every site 
belonging to the cluster. 

To facilitate the calculation of winding numbers, the HK representation is 
extended so that each pointer includes a distance, indicating the multiple of n 
from a given site to its representative. The path length from a given site to its 
root is then given by the sum over all distances along the path in the tree of 
representatives. Details about this can be found in |14| . 
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2.1 Identifying clusters 

When a torus is formed by the master from a rectangular lattice with open 
boundaries, periodic boundaries are applied across opposite sides of the rectan- 
gle. We call these boundaries the "seam" of the torus, which is shown in Fig. [I] 
The seam itself is a cross topology, so that it necessarily crosses all clusters 
which are not nomotopic to a point. Therefore, all winding clusters are encoun- 
tered during the gluing procedure and it is only there and then that one has to 
detect winding numbers. The seam is fastened like a zip that runs along pairs 
of sites. During this process clusters are detected to be winding and/or merge. 
If two sites either side of the seam are seen to be part of the same cluster prior 
to local gluing, because they are connected via another path (which may or may 
not cross the seam closed so far), a winding path has been detected, provided 
that the paths from the two sites to the same root differ. After the new path 
is normalized, two distinct cases need to be distinguished: 1) The cluster so 
far contains no winding path. In this case, the winding numbers are assigned 
to the cluster. 2) The cluster is already known to be winding, with winding 
numbers (a, b). If the new path (a', V) obeys (a', b') — (a,b), then the new 
path "runs parallel" to the old one. If the old and new paths differ, the cluster 
must have a cross topology: since both paths are normalized, they cannot be 
multiples of each other, because if there was a k such that (ka, kb) = (a' , b 1 ) or 
(ka 1 , kb') = (a, 6) one of the two paths would not be normalized. If (a, b) = (1, 0) 
then the new path must have b' ^ 0. Therefore (a, 6) intersects (a! ',&') thereby 
providing a "shortcut" and forming a path (0, 1). That is, (a' , b') taken forwards 
and (1, 0) taken a' times backwards gives (0, b'), which, normalized, is just (0, 1) 
(see below). 

The general case is that both paths wind in both directions. They are 
not multiples, so a/a' ^ b/b'. But then it is possible to find a path of the 
form (q, 0) with q =/= 0, namely by taking b' times path (a, b) and b times path 
(a' , b') backwards, which results in a path (b'a — ba', 0), with b'a — bo! ^ 0. All 
paths constructed by forward and backward movements may have their clusters 
deformed in such a way that the resulting path is composed of only forward 
movements. A path (q,0) with |q| > 1 must intersect itself, so that this path 
will always contain a path (1,0). 

The same construction provides a path (0, 1). Thus once again we arrive at 
the cross topology: if a cluster contains two incommensurable paths, then such 
a cluster has a cross topology. 

The recipe above provides a unique way to determine the winding numbers 
of a single cluster whenever different paths are identified in it. Similar principles 
apply to the entire configuration: The configuration itself is assigned winding 
numbers and a counter, keeping track of how often a particular type of winding 
cluster occurs. A case not mentioned above is that of two clusters, already 
known to be winding, that are merged across the seam. If the paths have the 
same normalized winding numbers, the counter is reduced by one. It cannot 
happen that they do not have the same winding numbers: It is topologically 
impossible that two clusters are both known to be winding with different winding 



7 



numbers before being merged at a later stage. Thus, if the first winding cluster 
is found, its winding number is assigned to the entire configuration and the 
counter is set to one. If a further, distinct cluster is found to be winding, it 
must wind with the same winding numbers, so that the counter is increased by 
one. If any cluster is found to have a cross topology, then all clusters will have 
a cross topology or be homotopic to a point. The cross topology is assigned to 
the entire configuration and the counter is set to one. 

Similar rules can be found for the Mobius strip. As explained above, a single 
cluster with winding number 1 can coexist with multiple clusters of winding 
number 2. Any new winding cluster found increases the total winding number 
of the configuration by its winding number. If two clusters merge, at least one 
of them must wind twice, reducing the total number of windings by 2. 

3 Results 

In this section we first list the key parameters of the simulation. Then we 
discuss the assumptions and estimates of the numerical errors associated with 
the study. In the following subsection we present the results for the torus which 
can be compared to Pinson's formulae Q and JHJ, asymptotes for which are 
presented in Section l3~4l In the following subsections, we present the results 
for the number distribution of multiple winding clusters and the results for the 
Mobius strip. 

In each simulation, we have determined Vn((cl, 6), n, r) as well as Pjv(0,r) 
and "P/v(X, r) for site and bond percolation at the critical point, i.e. occupation 
probability p (s) = 0.59274621 ^U] in site percolation and activation probability 
p lb) = 1/2 19 in bond percolation. 

Each simulation, parametrized by the system size N and the type of perco- 
lation (site or bond), consists of at least 10 6 realizations. A realization requires 
900 quadratic patches of size L 2 (L = 10, 100, 1000) produced by the slaves, 
which are joined together by the master, to form a lattice of size N = 900L 2 
and aspect ratio r, so that N = 30000 2 , 3000 2 , 300 2 . The resulting lattice is 
then "glued" together in different ways to form the different topologies. The 
14 different aspect ratios are: 30/30, 36/25, 45/20, 50/18, 60/15, 75/12, 90/10, 
100/9, 150/6, 180/5, 225/4, 300/3, 450/2 and 900/1. While a torus with aspect 
ratio r is topological^/ identical to a torus with aspect ratio 1/r, a Mobius strip 
can be glued along two different borders to form a band either of aspect ratio r 
or 1/r. Thus, 27 different aspect ratios are available for the Mobius strip. 

The same 900 patches are used for all different topologies and aspect ratios, 
but they are randomly permuted, rotated and mirrored before forming a new 
aspect ratio. The random number generator used is described in [20], and is 
especially suitable for parallel applications. 

The slaves, producing the 900 patches, themselves apply the appropriate 
boundary conditions, so that the simulation produces of the order of 10 9 samples 
for the torus and the Mobius strip for r = 1 and L = 10, 100, 1000. We restrict 
the presentation of these results to special cases. 
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3.1 Numerical errors 



In a numerical simulation an estimate p for the probability of the occurrence of 
a particular property (such as 2 winding clusters with winding number (1,3)) 
is measured as the average of an indicator function /(C) of the configuration C, 
which is 1, if the property is found in the configuration and otherwise. As 
all higher moments of / are p, the variance of p is then simply estimated as 
p — p 2 , so that the variance of the estimator of the probability is estimated as 
(p — p 2 )/(N — 1), where N is the number of (independent) realizations. 

It is worth mentioning a subtlety which can give rise to a systematic overesti- 
mation of the underestimation of the probability and vice versa: If the variance 
of the estimator for the probability p is estimated from the numerical data, and 
this data indicates a slightly smaller probability than the exact result, the vari- 
ance is likewise estimated slightly too small, if p < 1/2. Therefore, the relative 
numerical deviation of the numerically estimated probability from the exact re- 
sult will seem slightly too large. This becomes immediately clear if one wants 
to calculate the relative error for the probability for a configuration which actu- 
ally never occurred: the estimated probability and the estimated variance both 
vanish. Moreover, we note that the estimation of the variance itself acquires an 
increasing error as the probability approaches or 1. 

Vice versa, if the probability p < 1/2 is slightly overestimated from the 
numerical data, so too is the variance, and therefore the relative error is under- 
estimated. Correspondingly for p > 1/2, where overestimation leads to under- 
estimation of the variance and therefore to overestimation of the relative error 
and vice versa. However, since the aim of this work is to corroborate analytical 
work, all errors are based on the numerical results only, wherever possible. Only 
if the numerics give or 1 is the error based on the analytical result. If the 
analytical result is not available because of numerical convergence problems, no 
result is shown. 

As mentioned above, the same 900 patches are used multiple times to form 
rectangular lattices of different aspect ratios. Three different sets of boundary 
conditions (one for the torus, two for the Mobius strip) are then applied to the 
same rectangular lattice and the different measurements arc taken. Thus, there 
are correlations which, however, we do not explicitly account for, primarily be- 
cause we expect them to be very small. In fact, the correction to the error 
indicated depends on how many measurements are considered simultaneously: 
Assuming maximum correlations in n measurements, one could divide the sam- 
ple size by n, as if each quantity were based on its own subset of the sample. 
Accepting the averages calculated from the complete sample as good estimators 
of the averages calculated from each subsample, this procedure leads to a factor 
of \fn in front of each standard deviation. 

In the following we make use of the symmetry of the torus, 



V{{a,b),n, r) 
V(0,n,r) 
V(X,n,r) 



= V{{b,a),n,r 
= P(0,n,r- 1 ) 
= P(X,n,r- 1 ) 



-l 



) 
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i.e. the results of r and 1/r are not independent. 
3.2 Cross topology 

One surprising result in Pinson's paper |3j is that the probability of a cluster with 
cross topology (see the seam in Fig. 0), V(X,r), is identical to the probability 
of no winding cluster at all, P(0,r), i.e. the probability that all clusters are 
homotopic to a point. This is in perfect agreement with our numerical results: 
For N = 30000 2 Fig. H shows the deviation of V N (X,r) and Pjv(0,r) from the 
exact result (j^J for site and bond percolation*. 
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Figure 4: The deviation of ■pAr =3 oooo 2 (X, r) and 7 3 Ar = 3 ooo 2 (0, r) from V(X, r) (see 
(JHJ) in units of standard deviations for site (filled symbols) and bond (opaque 
symbols) percolation. 



It is obvious that there can only be one cluster with a cross topology, and 
that any other cluster must be homotopic to a point, i.e. if a configuration 
contains a cluster with a cross topology, there is no other non-trivial cluster. 

'There seems to be a subtle underestimation of V(0,r) and an overestimation of 'P(X, r), 
best seen in the slight segregation of triangles and circles. We are not sure yet what this effect 
could indicate. It is not present for N = 3000 2 and N = 300 2 . 
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3.3 Winding clusters 



Next we investigate the probability of at least one cluster with winding numbers 
(a, b), the exact expression for which was conjectured by Pinson to be QJ. Fig.[S] 
shows the relative deviation of the numerical results from the exact value for 
clusters with winding numbers (1, 0), (1, 1) and (1, 2) for N = 30000 2 and bond 
percolation. The reason why so many points seem to indicate no deviation at all 
is that the probability of certain types of winding clusters is extremely small. As 
a result, some rare types were not observed in the simulation and the resulting 
deviation from the exact result is approximately y/p/ y/N — 1, i.e. extremely 
small. Moreover, it should be noted that the numerical error of the evaluation 
of £[J increases as the probability approaches or 1. 
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Figure 5: The deviation of Pat = 3oooo 2 (( a > b), > l,r) from V((a,b),> l,r) (see 
0J) in units of standard deviations for (a, b) = (1, 0), (1,1), (1, 2) and site (filled 
symbols) and bond (opaque symbols) percolation. 



Fig.[S]shows a typical set of data for winding numbers (1, 2). It is clear that 
the probability for a (1,2) winding cannot be symmetric around r — 1. Above 
r = 2.0 it vanishes, because it becomes less and less likely for a cluster to wind 
around twice as the aspect ratio increases. Correspondingly, as r decreases it 
becomes less and less likely for a cluster to wind around once in the orthogonal 
direction. 
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All our numerical findings are fully consistent with JTJ for winding numbers 
(1,0), (1,±1), (1,±2) and (1,±3), in site and bond percolation. Other winding 
numbers occur too rarely to make any firm statements. 




Figure 6: The rescaled probability T^-Lsoooo 2 2), > l,r) for different aspect 
ratios. The rescaling factor is 10 4 . The straight line shows the analytical result 
©■ 

Percolation on a torus is mirror symmetrical, i.e. 

V((a, b), n, r) = V((a, -b),n, r) = V({-a, b), n, r) = V((-a, -b), n, r) . (7) 

While a cluster of type (a, b) is simply identical to a cluster of type (—a, —b), the 
windings (a, b) and (a, —b) are in fact distinguishable. Since these two quantities 
have the same probability, their comparison affords a consistency check upon 
the numerics, independent of any theoretical result or finite size corrections. We 
have passed this test successfully. 

3.4 Asymptotes 

To discuss numerical results for the asymptotic behavior of multiple winding 
clusters, we must first extract an appropriate functional form to fit against 
from (JTJ . A calculation closely related to the following was presented by Ziff for 
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the problem of crossing in two-dimensional percolation [211122) and subsequently 
used in [T4*] . 

It is convenient to slightly rewrite Eq. Q. Noting that all Z m ^ n enter ||TJ in 
the form Z a k.bk, it is reasonable to define 

Z((a, b);g,r) = J^J2 ex P (-^5 (~ + b 2 r)) , (8) 

which is just y^g/r #3(0, r), where r = ig(a 2 /r + b 2 r) and $3 is Jacobi's theta 
function |23| . By using transformations of the form ^ ( /22+1 = X); /( — Si /s; 
one has 

P((a,b),>l,r)= 1 (9) 
r;(e- 27rr ) 

x Qz(( 0> 6); 6, r) - Z((a, 6); 8/3, r) + ^((a, b); 2/3, r)J , 

given appropriate convergence of the sums in JSJ and QJ. For (a, 6) = (1,0) 
Eq. JSJ) transforms into a Riemann sum with mesh size 1/y/r, and for large r it 
converges to a Gaussian integral, i.e. lim r _ ! . o Z((a, b);g, r) — 1. That, however, 
leads to linv^oo exp(7rr/6) x in ©. This can easily be fixed using the Poisson 
summation formula (or the properties of $3 |23p. which gives 

lei igz 
and therefore for general (a, b) 

Z((a,b);g,r)= ] ^>xp (- * l 2 r) . (11) 
Va 2 + b 2 r 2 V g{a 2 + b 2 r 2 ) J 

From the definition of the Dedekind eta function 

00 

77(exp(-27rr)) = ]J (l - e - 2k7rr ) (12) 

fe=i 

one finds the following expansion for large r: 

V ~ 2 (exp(-2nr)) = e* r/6 (l + 2e~ 27rr + he~^ r + 10 e - 6jrr + . . . ) . (13) 
Using (|llfl for the large r expansion of @ at (a, b) = (1, 0) 

Z((l,0):g,r) = l + 2e- I T +2e- i2 T +2e- 92 T +2e- 162 f +... , (14) 
the task boils down to ordering terms: 
V((l, 0), > 1, r) = 1 - 2e - * ,rr + e - ^ + 2e~ 27rr - 4e"S 7rr + Ze^ r . . . (15) 
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n 


r range 


<7«((l,0),n) 


a((l,0),n) 


C< b >((l,0),n) 


a((l,0),n) 


1 


12/75.. 


.30/30 


0.811(3) 


1.568(3) 


0.810(3) 


1.566(3) 


2 


25/36.. 


.60/15 


0.856(3) 


7.771(12) 


0.868(4) 


7.813(14) 


3 


36/25.. 


.75/12 


1.364(10) 


19.02(4) 


1.368(11) 


19.02(5) 


4 


50/18.. 


.90/10 


2.05(2) 


33.91(8) 


2.10(2) 


34.05(10) 


5 


60/15.. 


.100/9 


3.76(6) 


53.94(16) 


3.71(7) 


53.82(19) 



Tabic 1: Multiple winding clusters (1,0). The numerical results 
7-V=30000 2 ((l,0),n, r) are fitted within the range indicated against lj2*0")l. 



This approximation has a relative deviation from the exact result of less than 
5 x 10~ 4 at r = 1 and less than 10~ 8 at r = 2. 

Using "P((l, 0), > l,r) ='P((0, 1),> 1, 1/r), the corresponding expansion for 
small r is now based directly on ©. Thus again for large r 

V((0, 1), > 1, r) = J (e-^ r - e~^ r - e~^ r + Ae~^ r . . .) . (16) 
V or V / 

At r = 1 the relative error of this approximation is better than 3 x 10~ 8 , which 
improves to about 10 -15 at r = 2. 
Similarly one finds for ^(X, r) 

Z o(f, 9i r) = 1 2 i((l, 0); / 2 3 , r) Z((l, 0); l/(/ 2 . 9 ), r) (17) 

which yields together with (|13l) and l|14|) 

p(X, r) = e - * 1 "" - e-i wr - + 2e~^ r - 2e^ 7rr . . . (18) 

again for large r. The relative deviation is less than 9 x 10~ 4 at r = 1 and about 
10~ 7 at r = 2. 

3.5 Multiple winding clusters 

Unfortunately, it is not straightforward to extend Cardy's qualitative arguments 
in the introduction of [7] for the existence of multiple spanning clusters. 

Multiple winding clusters with winding numbers other than (1,0) are very 
rare. In fact, we found 

^=10002 ((!. !). 2 > !) + P/tuxx* (C 1 ' 2 ' X ) = L40 ( 8 ) x 10_7 ( 19a ) 

^Liooo 2 ((1. !). 2 . !) + ^Liooo 2 ((1. "I). 2 > x ) = i- 50 ! 10 ) x 1(r? ( 19b ) 

based on the data produced at the slave nodes. Probabilities for higher multiples 
and other winding numbers are less than about 1 in 2 x 10 9 . 

From what has been said in Section l3~4l fsee also Eq. J2U) one might suspect 
that the probability of n distinct, simultaneously winding clusters with winding 
number (1,0) behaves in the limit of small r like 

P((1,0), n,r) « C((l,0),n)e- a « 1 '°)'™)/'"Vf . (20) 
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Figure 7: The rescaled probability 'P^L3oooo 2 ((^' 0)' n > r ) m tne form ln(7V(l — 
T 3 )) for different n. The dotted line shows the analytical result ijTjl. while the 
dashed lines are fitted according to H2U|) using the parameters shown in Tab. ^ 



The variable C((l,0),n) denotes the amplitude for this type of winding and 
a((l,0),n) is expected to be a second order polynomial in n [7j. Clearly, the 
additional factor y/r is "only" a logarithmic correction to the dominating expo- 
nential, but is in fact clearly visible in the numerics. 

Nothing can be derived from Section I3.4I concerning the large r limit, since 
the average number of winding clusters and the variance thereof become very 
large. 

Fig. H shows the estimated probabilities in the form ln(V/(l — V)) |2l I14|. 
These probabilities have been fitted against (|20[l in the small r limit, the results 
of which are shown in Tab. ^ Even though the choice of (1201 is somewhat 
arbitrary for n > 1, very good fits are obtained. These are shown by the dashed 
lines in Fig. [7| which represent the fitting results of Tab. ^fcd back into (|20|) . 
The main source of the numerical error is the fitting range, which is listed 
in the table. The range is bounded below by the smallest value of r that is 
supported by available data within reasonable error, and bounded above by the 
approximate value for r after which the asymptotic behavior terminates. The 
ambiguity of the fitting range is not reflected in the error bars, which indicate 
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only the statistical error. Therefore the exact result fitted to the function l|2U[l 
within the given interval should produce the values listed above. 

The dotted line in Fig.[7|shows V((l, 0), > 1, r), which asymptotically (small 
r) contains only a single winding cluster, i.e. 

Iim ffl,0),>l,r)-^((l > 0),l,r) =0 
r-o p((l,0),>l,r)+P((l,0),l,r) 

so that the amplitude C((1,0),1) and the exponent a((l,0),l) are known ex- 
actly from (|To|) : 

C((1,0),1) = V273 = 0.8164909 .. . (22a) 
a((l,0),l) = tt/2 = 1.5707963... (22b) 

which is in very good agreement with the results in Tab. ^ indicating that the 
fitting range chosen there is reasonable. 

Just as for wrapping clusters on the cylinder, one might be tempted to 
find a systematic dependence of C((l,0),n) and a((l,0),n) on n, such as an 
exponential and a second order polynomial, respectively. However, we were 
unable to identify these functions. Moreover, as the functional form l|20() already 
differs from the corresponding function for wrapping clusters on the cylinder, 
it is not surprising that no similarities were found between their exponents and 
amplitudes. 



3.6 Finite size effects 

The system sizes used are so large (9 x 10 8 sites) that one might be inclined to 
completely dismiss finite size corrections. Their study is nevertheless interesting 
for two main reasons. First, it is unknown a priori whether a system is actually 
large enough for finite size corrections to be neglected. Second, if they can be 
neglected, it is then interesting to investigate how strong the corrections are for 
smaller system sizes. 

The strength of finite size corrections has already been discussed in JS] for 
a simulation with the same parameters, for crossing, spanning and wrapping 
on the square lattice and the cylinder. There are, apart from Pinson's exact 
results, no estimates for probabilities specific to winding clusters on the torus 
that we know of. Fig.ODcannot be used as an indicator of high accuracy, because 
the aim of this paper is exactly to substantiate Pinson's analytical work and 
conformal invariance in percolation as such. 

However, it is very instructive to compare the numerical findings for differ- 
ent system sizes. Most remarkably, for the cross topology, the configuration 
homotopic to a point, as well as winding numbers (1,0), (1, 1) and (1,2), the 
distributions of the deviations from the exact results look essentially like Fig.0] 
and Fig. Eleven for N — 300 2 . The data corresponding to the latter are shown 
in Fig. EI 

We note that the deviations of the site percolation results from the exact 
values are not stronger than those of the bond percolation results. Thus we 
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to the corresponding data shown for N — 30000 2 in Fig. 03 



have also tested the estimate for p {s) = 0.59274621 |10|. which was derived on 
much smaller system sizes. 

3.7 Mobius strip 

The probability to obtain n windings on a Mobius strip is denoted in the fol- 
lowing by V{oo, n, r), where n encodes the number of winding clusters and their 
winding numbers as discussed above. Again, two different aspects can be in- 
vestigated: the probability "P(oo, > 1, r) and the probability "P(oo, a, r) for each 
individual a = 1, 2, ... . 

Fig. El shows the numerical results for bond percolation and N = 30000 2 , 
again in the reduced form ln(V/(l — T 3 )). The behavior is qualitatively similar 
to the one shown for multiple wrapping clusters on the torus, Fig- H Using the 
same ansatz as for wrapping clusters on the cylinder |14| . 

V{oo, n, r) w C(oo, n) e - a( - Q °^/ r , (23) 

the results shown in Tab. have been derived. The same precautions as in 
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Figure 9: The rescaled probability "PjvLsoooo 2 (°°> n ' r ) m ^ e f° rm m (?V(l — ^)) 
for different u. The dashed lines show the fitting results according to H23[l using 
the parameters shown in Tab. [21 The dotted line approximates the rescaled 
probability for r around 1 according to (|24bjl . 



Subsec. 13.31 apply to the ranges shown in the table. While a(oo, 1) can be 
conjectured to be 7r/3 with some confidence, the functional dependence of the 
exponents and amplitudes on n could not be found. Nevertheless, according to 
Fig-Elthe form l|23|) seems to work quite well. 

In the spirit of pQ we have tried to fit 'Pjv=30ooo 2 (°°7 ^ l> r ) against a third 
order polynomial in ln(r). The result 

V 1 ' JV=3000<M°°' - L > r )/ 
-0.6798(7) + 1.3525(9) ln(r) - 0.5648(9) ln(r) 2 + 0.2021(3) ln(r) 3 



In ( ^Uoooo 2 ^ 1 ^) ) « (24b) 
1 l-^Lsoooo^^ >!.»■) ' 



-0.6793(8) + 1.3545(10) ln(r) - 0.5689(11) ln(r) 2 + 0.2033(4) ln(r) 3 
is shown for bond percolation in Fig. as well. Similar to the results above, the 
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n 


r range 


C (s) (oo,n) 


a(oo, n) 


C (b) (oo,n) 


a(oo, n) 


1 


9/100.. 


.25/36 


0.982(3) 


1.0401(14) 


0.985(3) 


1.0433(17) 


2 


25/36.. 


.60/15 


0.671(2) 


5.832(9) 


0.676(3) 


5.849(11) 


3 


36/25.. 


.100/9 


0.780(2) 


15.62(2) 


0.774(3) 


15.54(3) 


4 


45/20.. 


.100/9 


1.081(11) 


30.93(10) 


1.099(13) 


31.08(12) 


5 


60/15.. 


.100/9 


1.44(6) 


50.6(4) 


1.51(6) 


51.1(4) 



Table 2: Winding clusters on the Mobius strip; winding n means n/2 winding 
clusters with winding number 2, if n is even, and (n — 1)/2 clusters with winding 
number 2 plus a single cluster with winding number 1, if n is odd. The numerical 
results 7 5 at = 3oooo 2 n i r ) are fitted within the range indicated against lj23|l . 

main source of error is not statistical, but systematic, namely in the choice of 
the specific function. Nevertheless, the numerical result (12411 will possibly serve 
as a reference for analytical findings. 

4 Conclusion 

Based on a large scale numerical simulation, this paper provides one of the first 
numerical confirmations of Pinson's analytical results for winding clusters on 
the torus, which are based on conformal field theory. It therefore also supports 
conformal invariance at the critical point. 

By rewriting Pinson's results, it was possible to derive some asymptotes that 
have hitherto only been derived for the flat topology |21U22| . These asymptotes 
have been used in the investigation of the probability of multiple, simultaneously 
wrapping clusters. 

A similar numerical analysis has been carried out for the Mobius strip, which 
still awaits analytical treatment. 
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